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Abstract 

The multivariate quantum g-Krawtchouk polynomials are shown to arise as matrix elements 
of “g-rotations” acting on the state vectors of many g-oscillators. The focus is put on the two- 
variable case. The algebraic interpretation is used to derive the main properties of the poly¬ 
nomials: orthogonality, duality, structure relations, difference equations and recurrence relations. 
The extension to an arbitrary number of variables is presented. 
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1 Introduction 

The purpose of this paper is to provide an algebraic model for the multi-variable g-Krawtchouk poly¬ 
nomials and to show how this model provides a cogent framework for the characterization of these 
orthogonal functions. The algebraic interpretation presented here is in terms of matrix elements of 
unitary “g-rotations” acting on g-oscillator states; these g-rotations are expressed as g-exponentials 
in the generators of several independent g-oscillator algebras. For illustration purposes, the focus 
will be put on the two-variable case. The algebraic model will lead to a natural and explicit derivation 
of the main properties satisfied by the bivariate quantum g-Krawtchouk polynomials: orthogonality 
relation, duality property, structure relations, g-difference equations and recurrence relations. How 
this approach generalizes directly to an arbitrary number of variables will also be explained. 

The standard (univariate) Krawtchouk polynomials form one of the simplest families of hyperge¬ 
ometric orthogonal polynomials of the Askey scheme [19], Ch. 9]. These polynomials of degree n in x, 
denoted by K n {x\p,N), have the expression 

K n (x-p,N) = (-m-N) n 2 F l 
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where 0 < p < 1 is a parameter and AT is a positive integeo In the above, p F q stands for the 
generalized hypergeometric series and 


(a) n = a(a + 1)• • • (a + n - 1), (a)o := 1, 


is the shifted factorial, or Pochhammer symbol Q] (1.1.2)]. The polynomials k„(x; p,N) are orthogonal 
with respect to the binomial distribution, i.e. 
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where r\ n ^ 0 are normalization coefficients. The Krawtchouk polynomials have a well-known alge¬ 
braic interpretation as matrix elements of unitary irreducible representations of so(3); in point of 
fact, most of their properties follow from that interpretation (see for example If20l|[28l Section 6.8]). 

In [26, 27] Tratnik presented a multi-variable extension of the Askey tableau, giving for each 
family of the scheme explicit expressions of the polynomials as well as of their orthogonality mea¬ 
sures; special cases of these multivariate polynomials occurred earlier, e.g. Bl8l . Much later in 
llllll . Geronimo and Iliev showed that Tratnik’s families of orthogonal polynomials are all bispectral: 
they exhibited the recurrence relations and the eigenvalue equations that these polynomials sat¬ 
isfy. Despite the fact that the multivariate polynomials defined by Tratnik are expressed in terms 
of one-variable polynomials, they are non-trivial and arise in a number of different contexts, see for 
example |9, 13, [18] EEL 22, 1241 . In two variables, the Krawtchouk polynomials defined by Tratnik 
have the expression 


K nu n 2 (xi,X2;pi,P2',N) = K ni (xi;pi,xi + X2)K n 2 (xi+ X2 - ni;p 2 ,N - ni), n\,U 2 = 0,1,.. .N, 


with n\ + n ,2 < N and are orthogonal with respect to the trinomial distribution 
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where { x ^ x ^ are the trinomial coefficients and where pi, p 2 are parameters. In IflOH . it was shown 
that Tratnik’s multivariate Krawtchouk polynomials arise in the matrix elements of the reducible 
unitary rotation group representations on oscillator states H- It was also shown that Tratnik’s 
Krawtchouk polynomials are in fact a special case of the multivariate Krawtchouk polynomials in¬ 
troduced by Griffiths in [12] and by Hoare and Rahman in |[15]. Similarly to the one-variable case, 
the properties of these multivariate functions follow from their group-theoretical interpretation IflOl : 
see also Ifl7l . 

Basic analogs of Tratnik’s multivariate orthogonal polynomials were introduced by Gasper and 
Rahman in 0. Like Tratnik, Gasper and Rahman explicitly defined their families of multivariate q- 
deformed polynomials as convoluted products of one-variable polynomials of the basic Askey scheme 
and proved their orthogonality relations. In Ill6ll . Iliev showed that these polynomials, which form 
a multivariate extension of the basic Askey tableau, are also bispectral. We note here that some of 

lr The polynomials K n (x\p,N) agree with K n (x\p,N) in 1191 Section 9.11] up to a normalization factor. In this paper, the 
uppercase K are reserved for the multivariate polynomials. 

2 Let us note that in ITOl and CD, a different definition of Tratnik’s two-variable Krawtchouk polynomials was used. 
The aforementioned algebraic interpretation is valid for both definitions. 
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these families of orthogonal polynomials, namely the multivariate g-Hahn and g-Racah polynomials, 
have already been seen to arise in algebraic and combinatorial contexts [23i [251. The two-variable 
g-Krawtchouk polynomials introduced in f5) have the expression 

(1) K niin2 (xi,X 2 ;a 1 ,a 2 ;N) = k ni (xi;a^ 2 ,xi + X 2 ',q) k n2 (xi + X 2 - ni',^ 2 ,N - n^q), 


with n\,n 2 are non-negative integers such that n\ + n 2 < N. Here k n (x,p,N;q) stands for the one- 
variable quantum g-Krawtchouk polynomials of degree n in q~ x which are defined as Ifl9l Ch. 14] 


(2) k n (x;p,N-q) = (-l) n (q N ;q) n q^ 201 
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where 0 < g < 1, p > g N and where N is a positive integei0. In the above 


(a;q) n = (1 - a)(l - aq) ■ ■ ■ (1 - aq n 1 ), (a;g) 0 :=l, 


stands for the g-Pochhammer symbol and r (/> s is the basic hypergeometric series defined as |4j 
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In this paper, we shall provide an algebraic interpretation of Gasper and Rahman’s multivariate 
(quantum) g-Krawtchouk polynomials and show how the properties of these polynomials can be ob¬ 
tained from that interpretation. On the one hand, this work can be considered as a multi-variable 
generalization of (8j, where the algebraic interpretation of the one-variable quantum g-Krawtchouk 
in terms of g-rotations was investigated. On the other hand, it can be viewed as a g-generalization of 
IllOll (in the Tratnik special case), where the interpretation of the multivariate Krawtchouk polyno¬ 
mials in terms of rotation group representations on oscillator states was studied. For simplicity the 
emphasis shall be put on the two-variable case, with the understanding that the approach directly 
extends to an arbitrary number of variables. 

The outline of the paper is as follows. In Section 2, the interpretation of the one-variable quantum 
g-Krawtchouk polynomials in terms of g-rotations, investigated in ®, is briefly reviewed. In Section 
3, the algebraic model for the two-variable quantum g-Krawtchouk polynomials is presented. The 
bivariate quantum g-Krawtchouk are seen to arise as matrix elements of unitary three-dimensional 
g-rotations acting on the state space of three mutually commuting copies of the g-oscillator algebra. 
The orthogonality weight and the duality relation for the polynomials are derived. In Section 4, 
new structure relations for two-variable g-Krawtchouk polynomials are found using the algebraic 
setting. In Section 5, the recurrence relations and the difference equations are obtained and are seen 
to coincide with those found by Iliev in Hl6l . The multivariate case is discussed in the conclusion. 


2 Review of the univariate case 

In this section, the algebraic interpretation of the one-variable quantum g-Krawtchouk polynomials 
as matrix elements of unitary g-rotations acting on the space of two g-oscillators is reviewed. Mild 
generalizations of the g-Baker-Campbell-Hausdorff formulas are also presented. 

3 The polynomials k n (x\p,N;q) agree with K% tm (q x ;p,N,q) of fl9l Section 14.14] up to a normalization factor. 
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2.1 The g-oscillator algebra and the Schwinger realization of U q (sl 2 ) 

Let 0 < q < 1 and consider two mutually commuting copies of the q -oscillator algebra with generators 
A+, A 0 and B 0 , B ± that satisfy the relations 

(3) [X 0 ,X ± ] = ±X ± , [X.,X + ] = q Xo , X-X + -qX + X- = 1, 

with X e {A,B} and where [X_, Y_] = 0 if X ^ Y. The algebra © has a standard infinite-dimensional 
representation on the orthonormal basis states 

I n A ,n B ) :=\n A )®\n B ), 


where n A and n B are non-negative integers. This representation is defined by the following action of 
the generators on the factors of the direct product states: 


(4) X+\ n x ) = 


'1 ~q 


nx +1 


1-q 


\nx + 1), X-\nx) = 


'i -q 


nx 


1 ~q 


■I nx ~ 1), 


■X’ol nx) = n x \n x ), 


where X = A or B. Note that in this representation, X + and X- are mutual adjoints and Xq is self- 
adjoint. It is easily seen that in the q t 1 limit, the representation defined by the actions © goes to 
the standard oscillator representation. 

The Schwinger realization of the quantum algebra U q (sl 2 ) is constructed from the two q -oscillator 
algebras © by taking El 


_Ao+Bo-l A 0 +E 0 -l Aq—Bq 

J+ = q 4 A+B-, J- = q 4 A-B+, Jo =--- 


Indeed, it is easily verified using the relations © that the operators J± and Jo satisfy the defining 
relations of I/ g (sl 2 ) which read 


Wo,J±]-±J ± , [J+,J-1 — 


q J o -q~Jo 

ql/2_ q -l/2- 


As can be verified directly using ©, the orthonormal q-oscillator states defined by 


(5) \n) N :=\n,N-n), n = 0,1,...,IV, 


where A is a non-negative integer, support the standard unitary (N + l)-dimensional irreducible 
representations of U q ( SI 2 ). 


2.2 (/-analogs of the BCH formulas and (/-rotation operators 

Consider the little and big q -exponential functions, respectively denoted by e q {z) and E q (z), which 
are defined as follows: 


( 6 ) 


j(z)= £ 


n t 0 {q ;q) n (z;q) 00 ’ 


\z\ < 1, 


00 q ( 2 ) z n 

and E q (z)=}^ - - — = (-z;q)^ 

n=0^T,q)n 


One clearly has e q (z)E q (-z ) = 1. In the following, two (/-analogs of the Baker-Campbell-Hausdorff 
relation shall be needed. The first relation is of the form 


00 1 n 

(7) E q (AX)Ye g (-Aq a X) = £ -- -[X,Y] n , 

re = 0 (<7;<7)„ 
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where [ X,Y\ n is defined recursively through 

[X,Y] 0 = 1, [X,Y] n+1 = q n X[X,Y] n -q a [X,Y] n X, n = 1,2,.... 

The second relation reads 

OO 1/1 

(8) e q (XX)YE q (-Xq a X) = £ -- ~[X,Y]' n , 

^o(q-,q)n 

where [X,Y~\' n is defined recursively by 

[X,Y]' 0 = 1, [X,Yl n+1 =X[X,Yl n -q n+a \X>Yl n X, n = 0,1,2,.... 

The identities (0, ® can be verified in a straightforward manner by expanding the q -exponentials 
in power series using ©. For a - 0, these relations are given in (3), among others. 

Let 0 be a real number such that \6\ < 1 and let Uab(9 ) be the operator defined as 


(9) UabW) = 

e)j 2 (d V Ao ) e q (0(1 - q)q- (Ao+Bo),2 A + B-} E q (-0(1 - g)g“ (Ao+Bo)/2 A_B + ) E^ 2 (-0 2 g“ s °). 

This operator will be called a “g-rotation” operator. This term, coined by Zhedanov in 112911 . comes from 
the fact that (7 ab( 0) can be viewed as a g-analog of an S'(7(2) element obtained via the exponential 
map from the algebra to the group; see [8} 29] for more details. As shown in 10, the operator ([9]> is 
unitary: it satisfies the relations (/'(/ = 1 and UU 1 = 1. The unitarity property of Uab(9) follows 
directly from the formulas 


e q (aA-B + )e q 


a,p 

,(1 ~q) 2 
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' a/3 
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where a,p,y,8 are constants or central elements. These formulas can be proved using the g-BCH 
relations 0 and ® as well as elementary properties of the g-exponential functions; see 0. 


2.3 Matrix elements and univariate quantum g-Krawtchouk polynomials 

Let us now recall how the one-variable quantum g-Krawtchouk polynomials arise in this setting. 
The matrix elements of the unitary g-rotation operator @ in the basis 0 are defined as 

(10) ^(G) = N (n\U A B(9)\x) N , 

where n and x take values in {0,1,..., N}. The matrix elements ( fTOl ) can be written as 0 

( 11 ) ^ x W) = (D^\d)cj^ ) {e)k n (x,B- 2 ,N-q), 

where k n (x,p,N;q) are the quantum g-Krawtchouk polynomials defined in 0. In ( fill ), the coefficient 
cr^ V) has the expression 


(12) <C'(0) = 


(-D re g ^ 


(q N \q)n \ 


g2nq-nN 

(9 2 q~ n -,q) n ’ 
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where 


M (q;q) N 

[n) q ( q-,q)n(q;q)N-n’ 

stands for the (/-binomial coefficient. The expression for a) x \0), which is defined as the matrix 
element n <0|C7 ab(9)\x) n, reads 


(13) oj ( x N) (d) := n (0\U ab (6)\x) n = (-1) 








{e*q- N -q) N 

,W 2 q- N -,q) x 
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From the algebraic interpretation flllD of the one-variable quantum q-Krawtchouk polynomials, one 
can derive the main properties of these polynomials (8J. In what follows, the duality relation 


(14) ®0) = .f 


(AO 

N—x,N—n 


W), 


which follows from the reality of the matrix elements ( IToT >. the identity U^ B (0) = U^jXO) and the 
observation that U^ B (0) = Uba(^), shall prove particularly useful. 


3 A model for the two-variable g-Krawtchouk polynomials 

In this section, the algebraic model for the two-variable q-Krawtchouk polynomials, which involves 
a three-dimensional (/-rotation acting on the states of three (/-oscillators, is constructed. The orthog¬ 
onality relation and the duality property of the polynomials are derived from the model. 

3.1 The model 

Consider the algebra generated by three mutually commuting q -oscillators with generators Xo, X±, 
X e { A,B,C }, satisfying the commutation relations (O and consider its representation defined by 
the actions 0 on the three-fold tensor product space spanned by the orthonormal basis vectors 
| riA,nB,nc) = \n A )®\nB)®\nc), where n\, ns and nq are non-negative integers. Let N be a positive 
integer and introduce the basis states \ni,n 2 )x defined as 

(15) \ni,n 2 )N :-\n 1 ,n 2 ,N-ni-n 2 ), n 1 ,n 2 e {0,1,...,1V}, n 1 + n 2 <N. 

For a given N, the states ( fl5l > support an (N + 1)(N + 2)/2-dimensional irreducible representation of 
the quantum algebra U q (sl 3 ) realized with three independent (/-oscillators; see fl4f Following ©, 
the unitary operator Uxy(S ) effecting the (/-rotation in the XY “plane” is defined as 

(16) Uxy(0) = 

e q 2 [0 2 q~ X °] e q (0(1 - q)g“ (Zo+Yo)/2 X+y_) E q (-0(1 - q)q“ (Xo+Yo)/2 Z_ Y + ) (-0 2 g“ y °). 

3.2 Matrix elements and bivariate q -Krawtchouk polynomials 

Consider the matrix elements 

(17) H^ ) „ 2 (xi,x 2 ;0,0) = N{ni,n 2 \U B c(<P) U AB (0)\xi,x 2 ) N , 

of the operator associated to a particular q -rotation in three dimensions. To ease the notation, 
E^ n2 (xi,X2;0,(/>) shall sometimes be written simply as E^ n2 (x\,x 2 ). An explicit expression for the 
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matrix elements ( 1171 ) can be obtained by inserting a resolution of the identity in between the two 
< 7 -rotation operators Ubc(<P) and Uab(8). This leads to 

N 

(18) Z™ n2 (x 1 , X2 ;9, ( p) = X N(nijTl^UBcitp)^,t)N N(s,t\U ab(6)\X1,X2)N■ 

s,t =0 
s+t<N 

Since the (/-rotation operator Ubc(<P) leaves invariant the first index when acting on the basis vector 
|rci,n 2 )iV, it follows from the definition dlOD and the orthonormality of the basis that 

(19) N(n 1 ,n 2 \U B c((P)\s,t)N = 6 nuS ^~ t ni \(p). 

Similarly, when acting on the state |xi,x 2 )jv, the operator Uab(S) leaves invariant the sum of the 
first two indices x\ + x 2 and one can thus write 

( 20 ) x(s,t\U ab(Q)\xi,X2)n = 3S + t,Xl+%2 ^S,X 1 (9)- 

Upon using d 19b and d20D in d 18b . the double summation cancels and one finds 

(21) H™ 2 (r 1) . 2 ;0^) = ^ 2) (0)Ci-« 1 W’ 

With the explicit expression dl lb of the one-dimensional (/-rotation matrix elements in terms of the 
univariate quantum g-Krawtchouk polynomials, one has 

(22) ^2 n2 (xi,x 2 ;6,(f)) = w i * 1 1+X2 \6)anl +X2) (6)k ni (xi,6~ 2 ,xi + x 2 ;q) 

x oj( ^+X 2 -n 1 ^° ( n ^ n ' \<P) kn 2 (xi + X 2 - n lt <p~ 2 ,N - m;q). 

Comparing d22l ) with Cl]), it is seen that the matrix elements H^ re2 (xi,x 2 ) can be expressed in terms 
of the two-variable q-Krawtchouk polynomials K„ 1 „, 2 (xi, v 2 ; ai, a 2 ;iV) with parameters cc\ = 9 and 
a 2 = <p. One must however verify that the matrix elements S^„ 2 (xi,x 2 ) factorize as a product of 
the matrix element n(0,0\Ubc( c I ) )Uab(8)\xi,X2)n, a normalization factor depending only on the de¬ 
gree indices n± and n 2 , and the polynomials K rei;n 2 (xi,x 2 ;ai,a 2 ;iV) themselves. A direct, though 
technical, calculation shows that this is indeed the case and that one can write the matrix elements 
“'?ii, 7 i 2 ^i’*^' 2 ) as 

(23) H^ n 2 (xi,x 2 ;0,0) = W^ 2 (0,0)Z^ re 2 (0,0)K ni;re 2 (xi,x 2 ;0,0;W), 
where W^ 2 ( 0 , 0 ) is given by 


(24) W^ X2 W,<p) = Ar(O,O|Usc(0)UAs(0)|xi,x 2 )iv 


= (-l) : 


,*2 


(_1)AT-X! ( g 2 -lya+x* q Q) 


< N 


*^1 5 %2 


(8 2 q\q)x 2 (<P z q\q)N- Xl -x 2 <P zr " q 


2 


2N-N(N+l)/2 


1/2 


where the normalization factor I.^ n2 (9,(p) has the expression 


(25) l<f n2 (0,0) = 9 -(V)-Cf)(-If 1 

(q\q)N 


(-D rai 


+n 2 


N 

ni,n2 


N(ni+n 2 )+2ni+n 2 —nin 2 —( n i)—{ n £) 


(qcp 2 -,q) n2 (q9 2 \q) [ 


/re i 


1/2 
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and where K„ 1 ) 7 j 2 (xi,X 2 ;ai,a 2 ;A 0 are the two-variable quantum g-Krawtchouk polynomials defined 
in l[T]l. In 0241 1 and 025D . stands for the q -trinomial coefficients 

'N\ ( q\q) N 

x,y) q ( q',q)x(q',q)y(q',Q)N-x- y ' 

Remark. Let us note here that in dl7D . the choice of the ordering of the two q -rotations to form a 
(^-rotation in three space is important. For example, consider the matrix elements 

N(n 1 ,n 2 \U AC (cl ) )U B c(e)\x 1 ,X2)N = 

which can be seen as q -analogs of the matrix elements leading to Geronimo and Iliev’s two-variable 
Krawtchouk polynomials II1111 . It can be verified directly that with this choice of ordering one has a 
factorization the form 

N(ni, n 2 \U AC (<l>)UBc(0)\xl,X2)N = n(0,0\UAc( c P'>Ubc(0)\ x 1, x 2)n Pn 1 ,n 2 ( x l,X 2 ), 
where P ni , ni (x with Po,o(u>* 2 ) = 1, are not polynomial in the proper variables. 


3.3 Orthogonality relations 

The orthogonality relation for the two-variable q -Krawtchouk polynomials immediately follows from 
the unitarity of the < 7 -rotation U B c((p)U ab(8), the orthonormality of the basis states and the reality 
of the matrix elements dl7D . One has 


(26) 


N(ni,n2\mi,m2)N = N(ni,n2\U B c( ( t , ')UAB(S)U' i AB (6)Ul c (cp)\mi,m2}N 

N 

= Y N(ni,n2\U B c((p)UAB(0)\xi,X2)N N(x 1 ,X2\U^ B (9)Ul c ((l))\m 1 ,m2)N 
* 1 >* 2=0 
Xl+*2<AT 

N 

~ ^m,n2^XljX2,0 ,(p) X2,0 ,(p) 

* 1 .* 2=0 

x\+x 2 <N 


Upon introducing d23l l in the above relation, one finds that the two-variable quantum q-Krawtchouk 
polynomials satisfy the orthogonality relation 


N 

(27) Y A ^ 2 ( 0 ,^)K/i 1; n2( x i,^2;0,</>;IV)K mi , m2 ( x i, x 2-,0,(p;N) = H ( ^n2(6,(p)8 nimi 8 n2 m2, 

xi,x 2 =0 

xi+x 2 <iV 


where A(xi,X 2 ) (iV) (0,0) = [W^ 2 (0 ,</>)] 2 and fT^ re 2 (0,0) = [X^ ra 2 (0,(/>)] 2 . From d26l ). it is easily seen 
that one can also write a dual orthogonality relation of the form 


(28) 


N 

Y h i ™n2( e ’^ K ni,n 2 (xi,X2]6,(p;N)K ni ' n2 (y 1 ,y2]0,(p;N) = H i x ^ X2 (6,(p)8 Xiyi 6 X2 y 2 , 

rei,re 2 =0 

rai+re 2 <iV 


where ( 0 , 0 ) = [!.)?’ ( 6 ,(p)] 


(N) 


and H™ 2 -' 


1 


-2 
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3.4 Duality 

The matrix elements dl7D satisfy a duality relation that allows to interchange the roles of the degree 
indices n.\ and with those of the variables indices xi and x<i- This duality relation is obtained by 
combining the duality relation dl4D and the expression d2ID for the matrix elements. One has 


(29) 


-(AO 


(xi ,x 2 ;0,cp) = &£ x 1 2 \0)Z 


(N- ni ) 

n2,xi+x 2 -ni 


(<£) = £ 


(xi+x 2 ) 
x 2 , xi+x 2 -n 1 


m 


(N- ni ) 

N-X1-X2, N-ni-n.2 


( (p ) 


AN-nA (th)f ( ' Xl+X2 ^ 

^N-xi-X2,N-n 1 -n 2 K V' ( ’x 2 ,Xi+X2-n 1 


( 0 ) 


-(AO 

“ N-Xl~X 2 ,X 2 


(N-ni-n 2 ,n. 2 ;(p,d), 


where we have used d2Tl ) in the first step and di ll ) in the second step; note that the parameters 6 and 
cp are interchanged in the duality transformation. In terms of the polynomials K raii „ 2 (xi,X 2 ;0,0;1V), 
the duality relation d29l ) becomes 


(30) (—l) X2 + ,l 2 qX\+x\+X\X2+Xi-N(Xi+X 2-1) (q® :Q)x2^ ( 1 ( P >Q')N-Xi-X 2 ^ 


\ (q',q) Xl 


K nun2 (x 1 ,X2;0,cl);N) 


_ n^+n^+nin2—ni—n2—N(ni+n2)\ >Q^n2 

— u 1 z - 

(q',q)N-ni~n 2 


^■N-xi-x 2 ,x 2 (N-ni -ri 2 ,n 2 -,(p,Q\N). 


4 Structure relations 

In this section, the g-BCH formulas ©, ©l and the algebraic model constructed in the previous sec¬ 
tion are used to obtain structure relations for the two-variable quantum g-Krawtchouk polynomials. 
Raising and lowering relations on the degrees and the variables are found. 


4.1 Structure relations with respect to the index n\ 

We first use the algebraic properties of the three-dimensional g-rotation operator Ubc(<I>)U ab(0) to 
obtain structures relations on the first degree index n\. 


4.1.1 Raising relation in n\ 

To get a raising relation with respect to the first degree index n±, consider the matrix element 
N-i(ni,ri2\A-UBc( ( p)UAB(0)\xi,X2)N ■ On the one hand, using Q) and ( fTTl ). one has 

1 1 _ qTIi +1 

(31) N-i(ni,n2\A-UBc(<l ) )UAB(0)\xi,X2)N = y — T^q — 

On the other hand, one can write 

(32) N -i{n 1 ,n2\A-U BC {(p)U AB mxi,X2)N = N-i{ni,n2\U B ci ( P)U AB iq 1/2 0) A_|xi,X2>Ar, 
where A_ is of the form 

(33) A- = u\ B {q~ ll2 d)Ul c {(p)A - UbcWUabW) = u\ B {q~ ll2 d) A_ U AB (8). 

A direct calculation using the expression d 16D and the formulas (0 and I© yields 

(34) U\ B (q~ 1/2 e) A_ U AB (9) = \J 1 - e 2 q- 1 q~ B ‘>A _ + 0g“ 1/2 g (A °“ s ° )/2 S_. 
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To obtain d34l ) from the g-Baker-Campbell-Hausdorff relations, one should recall that U' XY = ^XY' 
Upon substituting ( l34l ) in ( f32l ). using the actions ft]) and comparing with ( HTTb . one finds 


(35) yJl-q n i + 1 E { ^ ln 2 (x 1 ,x 2 ;9,(f))=yJl-e 2 q ** 1 y/l-q x iE^ n 1 J(x 1 -l,x 2 ;q 1 ! 2 d,(p) 

+ 0g tel “ X2)/2 V *x-V{ Xl ,xz - 1 ;<T 1/2 e,0). 

Using the expression d23l ) of the matrix elements in terms of the two-variable g-Krawtchouk polyno¬ 
mials, the relation d35D takes the form 

, 0 „, q- ni K ni+ltn 2 (x 1 ,X 2 ;q y 2 6 ,(P;N) = q- (x ' +x *\l-q x ')a-0- 2 q x *)K nun 2 (x 1 -l, X 2 ;0,(P;N-l) 

(36) 

- (1 - g 2 )K„ 1>re2 (xi,x 2 -l; 0 ,(/);A(-l). 

4.1.2 Lowering relation in«i 

A lowering relation on the first degree index n \ can be obtained by considering instead the matrix 
element N+i(ni,n 2 \A+UBc(<P)UAB(q~ 1 / 2 9 )\xi,X 2 )N■ One has 

(37) JV+l(«1.^2|A + Usc(0)UAs(g _1/2 0)|Xl,X2>iV = 1 - q ^ni-l,nS Xl A2',q~ ll2 9,<p), 

and also 

(38) N+i( n i,n2\A + UBc(<l ) )UAB(Q- 1/2 9)\xi,X2)N = N+i(ni,ri2\UBc( c P')UAB(9')A + \xi,X2)N> 
with 

(39) A+ - u\ B W)Ul c {<p)A + U B c(<p)UAB{q- ll2 0). 

It is directly seen that d39l) is the complex conjugate of d33l) . Consequently, we have 

(40) A+ = A + \J 1 - 6 2 q- l q~ B ° + 0g“ 1/2 B+g (A °“ So)/2 . 

Upon substituting d40l ) in d38l ). using the actions Q and comparing with < T37T >. one finds 

(41) v / W rr S^ ln2 (x 1 ,x 2 ■,0,<p)=yjl-q*i + 1 y/l-6 2 q-x* sj^ + l ,x 2 ;g 1/2 0,^) 

+ 0g (xi_X2)/2 y / l. - q x 2+1 "^(xi.xa + l-,q 1 / 2 d,(p), 

where we have taken d —>• q 1/2 d. Using ( l23l ). this gives 

Q n 1(1 — O ni ) 

(42) - — -K rei _i ;n2 (xi,x 2 ;0,^;A/ r ) = g* 1+2 K„ 1>n2 (xi,x 2 + l;q 1 / 2 0,(p;N +1) 

-g" 1+2 K ni> „ 2 (x 1 + l,x 2 ;g 1/2 0,</);Ar+ 1). 

4.2 Structure relations with respect to the index 

We now look for structure relations on the second degree index 
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4.2.1 Raising relation in ri 2 


To get a raising relation with respect the second degree index ri 2 , consider the matrix element 
N -i(n 1 ,n 2 \q Aol2 B-U BC ((l))UAB(Q)\xi,X 2 )N- One has 


/ 1 _ q 112 + 1 

(43) N-i(n 1 ,n 2 \q Ao,2 B-U B c(cp)U AB (d)\x 1 ,X 2 )N = <? rei/ V 1 _ 2 +i( x i A 2 ;d,</>)- 

Proceeding as in the previous subsection, one writes 

(44) N -i(ni,n2\q Ao/2 B-U B c(^)UAB(0)\x 1 ,X2)N = N-i{ni,n 2 \UBc(q~ V<2j( p)U ab(S)B-\x\,X 2 )n, 
where 


(45) B- = U\ B {6)Ul c {q- in <p) q~ A ° /2 B _ Ubc^UabW). 

Upon using d34l ) with the As replaced by Bs and the Bs replaced by Cs, one finds 

(A 0 +B 0 )/2 


1 -0 2 g- 1 g- Co g _Bo/2 fi-+0g _iy2 g _Co/2 C- \U AB m. 




It is obvious that q A0+B ° commutes with U AB (S), and similarly for all the terms involving the gener¬ 
ators Cq, C+. Using the conjugation formula 


U\ „(9) q~ B oI2 B _ U AB (9) = q~ B °' 2 B _ J 1 - 6 2 q ~ B ° - d q ~ B ° A_q 


-A 0 /2 


J AB 

which can be obtained by a straightforward calculation with the help of the g-BCH formulas (7J), ©, 
one obtains 


(46) B. = q 


_ <A 0 +B 0 )/2 


1 ~(p 2 q 1 q c ° 


, q~ B o ,2 B_ J 1 - e 2 q~Bo - 9 q- Bo A-q~ Ao/2 ^ 


+ <Pq- ll2 q- C ° l2 C _ 


Upon inserting A46D in H44D . and comparing with ( 143D , one finds 


(47) g ni/ V 1 - q^~.™ n2+l { Xl ,X 2 -,e,cp) = <Pq x ^q~ N/2 ^ 1 - q N ~* H^W 2 ;0,g- 1/2 0) 


q- 1/2 dq- X2,2 Jl -(p 2 q x i+x 2 -N-l ^/i _ qXl eW-»( Xi - 1 ,x 2 ;9,q- 1/2 (p) 


+ q Xll2 sJl-q x *\J 1 - 0 V * 2 yj 1 - ^qM-N-l - 1 ;d,q 1/2 (f>). 

Using the expression ( l23l ). one obtains 

(48) q N ~ n2 K nun2+1 (x 1 ,x 2 ;6,q 1/2 <P;N) = (1 - q x ')( 1 - (TV"* 1 ”* 2 ) K nij „ 2 (xi - 1 ,* 2 ;0,0;1V- 1) 

+ ^(1 - q X2 )(l-(p- 2 q N - Xl - X2 )K nun2 (x 1 ,X2 - l;d,<f>;N-l) 

+ q Xl+X2 { 1 - g*"* 1 -* 8 ) K niiB2 (*i,* 2 ;0,0;N-1). 


The relations d36l ) and d48D can be combined to generate K„ 1)re2 Ofi > ^ 2 ; «i, « 2 ; N) recursively. 
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4.2.2 Lowering relation in ri 2 


In order to obtain a lowering relation in the second degree index n 2 , one considers the matrix ele¬ 
ments N+i(ni,Ti 2 \B + q Ao/2 UBc(q- 1/2( P')UAB(Q)\xi,X 2 )N- One has 


(49) N+i(ni,n 2 \B + q A ° / 2 U B c(Q uz <p)U ab { 8 )\x\,X 2 )n = <7 


- 1 / 2 , 


rei/2 


'1 ~Q ni 


1 _ q ~n,]n 2 -l( x l, x 2-,0,q ll2 <p\N), 


and also 

(50) N+ 1 (ni,n 2 \B + q Ao, 2 UBc(q~ 1 / 2 ^)U AB ( 6 )\xi,X 2 ) N = N+i(ni,ti 2 \U B c((p)U AB ( 8 ) B+\xi,X 2 )n, 
where B + is given by 


(51) B + = U\ B {6)U BC {<p)B + q^ U BC {q-^<p)U AB {8). 

It is seen that 05ID is the complex conjugate of 045D . Taking the complex conjugate of ( 146 D . one finds 

-C 0 / 2 I (A 0 +S 0 )/2 


7 A 0 /2 ; 


.-1/2, 


(52) B + = 


1 - d 2 q~ B °B + q~ B ° /2 - dq- Ao/2 A + q- Bo J 1 - <p 2 q- l q ~ c ° + (Pq~ V2 C + q 


Proceeding in the same manner as before leads to the relation 

(53) g^y/T^ S^ ||2 _ 1 (*i > * 2 ;0, 9 - iy2 ^) = 

-Oq~ 1,2 q~ X2l2 \J 1 - cp 2 q~ 1 q x i+ x i~ N yj 1 - q*> 1+ 1 + l,x 2 ;0,(p) 

+ q Xli2 \J 1 - <p 2 q- 1 q Xl+X2 ~ N \J 1 - 8 2 q~ 1 q~ X2 'Jl - q X2+1 + 1 ; 0 , 0 ) 

+ <pq~ ll2 q Xl+X2 q ~ N/2 \Jl- q N ~ x 1 -^+ 1 ,x 2 ;0,(p). 


In terms of the bivariate g-Krawtchouk polynomials, this relation amounts to 

rei+re 2 Q_„re2) 

(54) - - -K n n2 - 1 {xi,X2]9,q~ m <p\N) = -6 2 q Xl K re „ 2 (xi + l,X 2 \d ,(p\N + 1) 

</) z 

+ e 2 q x m - e- 2 q X 2 + 1 )K ni , n 2 ( Xl ,X2 + 1 -,e,<p-N+ l) + ? X1+X2+1 K ni!n2 (x 1 ,x 2 ;8,(P;N + 1 ). 

4.3 Structure relations with respect to the index x\ 

Proceeding in a “dual” fashion, one can obtain structure relations in the first variable index x\. 


4.3.1 Lowering relation in x\ 

In order to find a lowering relation in first variable index x\, one considers the matrix element 
n( ni, i%2\U B c((p)U ab ( 6 ) q~ A ° / 2 A-\xi,X2)N+i- One the one hand, one has 

(55) iv(n 1 ,n2lt/Bc(0)^AB(0)g“ Ao/2 A4x 1 ,x 2 ) iV+ i = 9 - (xi - 1)/2 ^/^^H^ 2 (xi-l,x 2 ;0,0;lV). 

On the other hand, one can write 

(56) N( n i,n 2 \U B c((p)U AB ( 8 ) q~ A ° / 2 A-\x\,X 2 )N+i = N(n\,ii 2 \ A_ U B c((p)U ab { 8 )\x\,X 2 )n+i, 
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where 


(57) A- = UbcWJabW) q~ A ° /2 A _ U^mU^). 

To compute ( f57l ). one first uses the conjugation formula 


(58) U A B(0)q~ Ao/2 A-U\ b (B) = q~ A ° 12 A_ \Jl-9 2 q~ A °- 9 q- A °B_q- B ° /2 , 

which is obtained straightforwardly from the g-BCH formulas Q) and I®; this lead to 

(59) A- = U bc (<I>) 


q~ Ao,2 A_ yjl -9 2 q- A °-9q- 1,2 q- Ao q- Bo/2 B- 
Using < f58l > again with A replaced by B and B replaced by C, one finds 
(60) A- = \q~ Aol2 A. J 1 - 9 2 q ~ A » - 9q~ 1/2 q~ A ° "~ BJ2 : 


u' BC w- 


B- v l-(p 2 q 


2 n ~Bn 


(pq 


-1/2 -Bo -C 0 /2 


C- 


Upon substituting d60l ) in ( l56l ). using the actions 0 and comparing with d55l ). one finds 

(61) q- (x i- 1V2 ^T^E™ n2 (x i-l,* 2 ;0,0) = 

q~ n 1/2 \/l - g rei+1 \/1 - B 2 q- 1 q~ n ^ x 2 -,8,(p) 

- 0<r 1/ V rei <r re2/ V 1 - q n * +1 ^ 1 - W" 2 ’ x 2',B,(p) 

+ 00g-V (ni+B2)/ V iW2 

In terms of the polynomials, this amounts to 

(62) -^-K rei , re2 (xi - 1 ,x 2 ;9,(f);N) = q~ 2ni K ni+ i y n 2 (xi,x 2 -,6,(p-,N +1) 

+ cp 2 q~ ni ~ 2n2 ~ 1 K nitn2+1 (x 1 ,x 2 -,9,(f);N + l)-(p 2 q~ N ~ 2 {l - q N ~ n ^ r 


(p). 


— 77,0 + 1 \ 




4.3.2 Raising relation in x\ 

A raising relation in the first variable index x\ can be obtained by considering the matrix element 
N(ni,n 2 \UBc(<p)U ab( 8) A+q- Aol2 \xi,x 2 )N-i- One has 

(63) N (ni,n 2 \UBc((p)UAB(9)A+q~ Ao/2 \xi,x 2 )N-i = q~ Xl/2 ]J — E ( ^ n2 (xi + l,x 2 ;d,(p), 
and also 

(64) n (ni,n 2 \UBc(<P)U ab(B)A + q A ° / 2 \x\,x 2 )]s[-i = N(n\,n 2 \A + TJbcWUab(Q)\x\,x 2 )n+\, 
where 

(65) A+ - UbcWUabW) A + q- Ao,2 u\ B {9)Ul c (cp). 
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The expression for d65l ) is obtained directly by taking the complex conjugate of ( 1601 ). Substituting the 
result in A64D and comparing with ( 163D . one finds 


(66) <T* i/2 /l - + 1 ,x 2 ;d,cp) = q m q- n ^yj 1 -0 2 g-»i y/l^ 

-dq-^q-^y/T^y/l-ilPq-"* E^^Xi MM) 

+ e<Pq- 1/2 q- (N+ni+n * )/2 yJ 1 - qN-m-n 2 

Using the expression d23D for the matrix elements, this relations translates to 

(67) K„ 1>re2 (xi +1 ,X 2 \ 6 ,(p\N) = g _1 (l - g rei )(l -0“ 2 g rei )K„ 1 _i ire2 (xi,x 2 ;0,0; N- 1) 

+ q- 1 - ni a-q n2 )a-(p~ 2 q n2 )Kn 1 ,n 2 -l( x l, x 2-,e,(p;N-l) + q- ni - n2 K nu n 2 (x 1 ,X2-,e,(P;N-l). 

4.4 Structure relations with respect to the index 

We now obtain structure relations with respect to the second variable index X 2 - 


4.4.1 Lowering relation in X 2 

In order to obtain a lowering relation with respect to the second variable index X 2 , consider the 
matrix element at (n\,ri2 \Ubc^Uab(6) Q C ° /2 B~\xi,X 2 )n+i • One has 


(68) N (n 1 ,n 2 \U BC (cl))U AB (e) q Co, 2 B.\ Xl ,x 2 ) N+ i = gW-*i-**+iV2 


ct( ao 


l-g 




. 1/2 c 


and also 

(69) N(ni,Tl 2 \UBc( ( P)UA b( 6 ) q C ° i ‘ 2 ‘B-\Xi,X 2 )N+l = N(n\,ri 2 \B- Usc(</>)UAB(g 1/Z 0)|xi,X2}jV+l, 
where 

(70) B _ = UbcWUabW) q Col 2 B- U\ B (q in e)Ul c (cp). 

To determine the expression for B-, one needs the conjugation relation 


(71) U AB (0)B- U^ AB (q 1/2 d) = y l-0 2 q~ A °B- + 0g (So-Ao)/2 A_, 

which is again obtained from the g-BCH relations. Using this relation, d70D becomes 

(72) B- = UbcW {g (So+Co)/2 [y/ 1 - 9 2 q~ A ° q~ B ° l2 B- + 0g“ Afl/2 A_] | t/J c (0). 

Upon using d58l ) in the above (with A-+B and B —>■ C), one finds 

2 }{\/l ~9 2 q 


(73) B- = \ q 


JB 0 +Co)/2 | 


y-Ao 


,~B 0 /2- 


B-Jl-(p 2 q B °-(pq B °C-q 


-Co/2 


+ 9q 


-A 0 /2 


A_ . 


Substituting the above expression in d69l ) and comparing with d68l ). there comes 

(74) g w -*i-*2 + D/2 v /rr^ sJJ #2 (x 1i x 2 - 1 = 9q N/ 2 q~ n ^l - g" 1+1 E™ a (xi,x 2 ; g 1/2 0,</>) 


+ q 


(N-n\-ni)l 2 


1 -9 2 q- n ^l-q n * +1 ^l-<p 2 q- l q- n * H^^(xi.xajg^M) 


-(Pq- 1 / 2 q- n 2 / 2 Jl-9 2 q-^Jl-q N ~ n ^ +1 E^t 1 J(x 1 ,x 2 ;q 1 / 2 e,(f)). 
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In terms of the polynomials K ni>n2 (xi,x 2 ;9,(p;N) this relation is of the form 


(75) 


Q l-Xl~X2 

—- (1 - 9 X2 )K ni> „ 2 (x i,x 2 - l;B,<p;N) = g 1-2 ” 1 K ni+ ^ n2 {x 1 ,x 2 -,q m B,(p\N + 1) 

+ <p 2 q~ ni ~ 2n2 ( 1 - 9~ 2 q ni )K ni! n 2 +l(xi,x 2 ;q 1/2 0 ,(p;N + 1) 

- <p 2 q- N ~\l - g tf-»i-» a +lxi - 9- 2 q ni )K ni ^{x l ,x 2 -q ll2 9,<p-N+ 1). 


4.4.2 Raising relation in x 2 

To obtain a raising relation in the second variable index x 2 , one considers the matrix element 
N(ni,n 2 \UBc((p)UAB(q 1/2 9)B + q c ° /2 \xi,x 2 )N-i- On the one hand, one has 

(76) Ar(«i,« 2 |(/Bc(0)(7As(g 1/2 0)5 + g Co/2 |xi,x 2 > A r-i = g w - Xl - X2 - 1)/2 W^^- H( ^2^i^2 + l;0^), 

and on the other hand one can write 

(77) N( n i,n 2 \UBc( c l ) ')UAB(Q ll2 9)B + q Co/2 \xi,x 2 )N-i = n(tii, n 2 \B+ Ubc(<P)Uab(9)\xi,x 2 )n-i, 
where 

(78) B + = U B c((p)U A B(q 1,2 9)B + q Co/2 U AB (6)U BC (<p). 

Since B + =B 1, one directly finds 


(79) q W-xi-x»- ivy 1 - gx 2+ l Z™ n2 ( Xl ,x 2 + l-,q 1/2 9,(p) = Bq^q-^q^^l^E^Jx^B^) 
Y (JV-ra 1 -/i 2)/2 / 1 _ ^ q - n2 ./i _ 82 q - ni sW-U( Xl ,x 2 ;d,(l>) 


+ q 


( pq-^q-^^ 1 _ q N-n 1 -n 2 ^ 1 _ 92 q-n 1 eW-1>(*i,* 2 J0.0). 


In terms of the polynomials K rei>n2 (xi,x 2 ;0,</>;lV), this gives 

(80) K„ 1>re2 (xi,x 2 + l;g 1/2 0,<();iV) = g _1 (l-g ni )K„ 1 _i ;re2 (xi,x 2 ;0,</);iV-l) 

+ g- rei - 1 (l-0- 2 g' l2 )(l-g re2 )K„ 1;re2 _ 1 (xi,x 2 ;e,0;iV-l) + ( ? - (rei+re2) K nijre2 (xi ) x 2 ;e,0;iV-l). 


5 Bispectrality 

In this section, the algebraic model for the two-variable g-Krawtchouk polynomials is used to derive 
their recurrence relations and their eigenvalue equations. These relations coincide with those of 11611 . 

5.1 Difference equations 

To derive the two eigenvalue equations that determine the polynomials K ni> „ 2 (xi,x 2 ;ai,a 2 ;.ZV), one 
could find suitable combinations of the raising and lowering relations on the degree indices n i and 
n 2 obtained in the previous section. However, it will prove more instructive to proceed directly by 
finding conjugation formulas for the appropriate operators. 
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5.1.1 First difference equation 

Consider the matrix element N(n\,ri 2 \q Ao Ubc(<P)Uab(8)\xi,x%)n- On the one hand, one has 

(81) N<n 1 ,n2\q A ° U B c(<P)Uab(Q)\x\,X2)n = q ni ^n 1 ]nS x ^X2\d,(p). 

On the other hand, one can write 

(82) N(n\,n2\q A ° Ubc(<P)Uab(9)\xi,X2)n = N(ni,n2\UBc(<P)U ab( 8) q A °\xi,X2)N, 
where 

(83) ^ = U\ B (9)Ul c {<p) q A ° UbcWUabW) = U\ B (0) q A ° U AB (0). 

The conjugation U^ B (d) q A ° Uab(9 ) can be performed with the help of the g-BCH formulas. As this 
computation is not straightforward, we present the details of the calculation. Using the explicit 
expression ( fl6l l for the q -rotation operator, one writes 

q^° = e^ 2 (0 V- 80 ) e q (0(1 - g)g“ (Ao+So)/2 A_B+) E q (-0(1 - q)q- (A ° +B ° )l2 A + B-} 
x E]j 2 (-0 2 g _A °j |c ? Ao }e^ /2 (0 2 c ? “ A<) ) 

x e q (0(1 - q)q- (Ao+Bo),2 A + B- ) j E q (-0(1 - g)g“ (Ao+So)/2 A_B + ) E^ 2 (-0 V 8 °) . 
The first conjugation is trivial and one can immediately write 

q^> = e]' 2 (0 V 80 ) e q (0(1 - g)g“ (Ao+8o)/2 A_B+) 
x E q (-0(1 - g)g“ (Ao+8o)/2 A + B_) q A ° e q (0(1 - g)g“ (Ao+8o)/2 A + B_) 

x E q (- 0(1 - g)g“ (Ao+8oV2 A_B + ) E^ 2 (-0 V 8 °). 
The next step is performed using the identity q A ° e q (XA+) = e q (qXA + )q A °; one then finds 
q^ =e 1 c { 2 (0 2 g- 8 °) J 1 J2E 1 q /2 (-e 2 q- Bo ), 
where 


(84) J± = e q (0(1 - g)g _(Ao+8o)/2 A_B + j l-0(l-g)g- (Ao+8o)/2 A+B_ E q (-0(1 - q)q- (Ao+Bo)/2 A.B, 

and 


(85) J 2 = e q [e(l-q)q- {Ao+Bo),2 A-B- 


E q (-0(1 - q)q- (Ao+Bo),2 A.B + ) j. 


The expression for A84D can be obtained from the conjugation formula 

a A _ a 1 


e g (aA-B + )A + B-E„(-aA-B + ) = A+B- + ■ 


9 ' + + ' ?V ’ (1 ~q) 2 

which follows straightforwardly from ([8]). One thus has 


q A ° ~ 


-.Bn 


(l-q) 2 1-aA-B, 


Ji = l-9(l-q)q 


-(A 0 +B 0 )/2 


A+B- + 


eq 


-(A 0 +Bo)/2 


-.Ao 


(1-g) 


Q q-(Ao+Bo)/2 


(l-q) 1-6(1-q)q~^ +B ^ 2 A-Bo 


-.Bo 
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One easily finds 


J 2 = q A ° [l-0(l -q)q~ (Ao+Bo)/2 A.B + 

Performing simplifications on the product of J\ and J 2 using the commutation relations one finds 

> - e* /2 (0 2 g“ B °) |g A °[ 1 - 0 2 g“ B °] -6(1- g)g (A °“ Bo)/2 [l - 0 2 g“ s °]A_ J B + 

- 0(1 - q)A + B_q (A °- Bo)/2 + q A ° (d 2 - g _1 0 2 q~ B ° + 0 2 g“ 1 g“ (Afl+Bo) ) j E]j 2 (-0 2 g“ B °). 
The last conjugation is easily carried out; the final result is 
(86) q^> = q A o [ i _ (1 + g- 1 )0 2 g“ B o +e 2 + q- 1 6 2 q- {Ao+Bo) 

- 0(1 - g)g (Ao “ Bo)/2 \Jl-9 2 q~ B °A-B + - 0(1 - g )A+5 _ \J 1 - 0 2 q~ B ° g (A ° “ Bo)/2 . 

Upon using A86H in d82D . using the actions Q and comparing with d81D . one finds that the matrix 
elements E^ n 2 (xi,x 2 = E^ n 2 (xi,x 2 ) satisfy the eigenvalue equation 


(87) q ni E (Af 1 > n 2 (xi,x 2 ) = [g Xl (l + 0 2 ) + q~ 1 6 2 q ~ X2 - (1 + g _ 1 )0 2 g Xl_X2 ] E™, a (*i,* 2 ) 


— In2 —X 2 


-l\n2 Xi ~x 2 ] 77 CAT) 


- 0 g _ 1 g (xi_X2V2 


(1 - g X2+1 )(l - g Xl )(l - 0 z g 


l 2 -(^2 + 1 ) 


1/2 


2 + 1) 


- 0 g (xi “ X2)/2 [(1 - 0 2 g _X2 )(l - g X2 )(l - g Xl+1 )] 1/2 (*i + l,x 2 - 1 ). 


Upon using the expression d23D for the matrix elements, one finds that the two-variable quantum 
g-Krawtchouk polynomials K ni) ,j 2 (xi,:r 2 ;0,</>;AO satisfy the 3-point eigenvalue equation 

( 88 ) q ni K nun 2 ( Xl ,x 2 ) = [g Xl ( 1 + 0 2 ) + g- 1 0 2 g “ X2 - (1 + g- 1 ) 0 2 g Xl “ X2 ] K„ 1 |B 2 (*i,* 2 ) 

+ (1 - 0 2 g“ X2 “ 1 )(l - g Xl ) K, ll>re2 (xi - 1 ,X2 + U + 0 2 g* 1_ **(l - g X2 ) K ni> „ 2 (*i + l,x 2 -1). 


5.1.2 Second difference equation 

Another difference equation for the matrix elements Z^ n 2 (xi,x 2 ) and the two-variable Krawtchouk 
polynomials can be obtained from the matrix element N(ni,n 2 \q Ao+B ° Ubc(<P)U ab(9)\x\,x 2 ) at ■ One 
can write this matrix elements as 

(89) jV<rai,rc2lg Ao+S ° Ubc(<P)U Ab(8)\X\,X2) N = g rei+re2 E^ ) n2 (:r l> x 2',0,<p)- 
but also as 

(90) N( n i> n 2 \q Ao+B ° Ubc(<P)U ab(6)\xi,x 2 )n = N(ni,n 2 \UBc(^)UAB( 6 ) q Ao+B °\x i,x 2 )n, 
where 

(91) gCT'o = u\ B (6)Ul c (<p) g Ao+B ° Ubc((P)Uab(9). 

Using the result dMt with the As replaced by B s and the B s replaced by C s, one has 


(92) qAo+B 0 — U^ AB (d) g Ao+s o | l-(l + g- 1 ) 0 2 g- C ''+ 0 2 + g- 1 0 2 g-( B o+Co) 

- (p(l - q)q- (Bo+Co)l2 J 1 - (p 2 q~ c °B-C+ - </> g _1 (1 - g)B+C_ J 1 - 0 2 g- c ogABo+C o )/2 I 
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To complete this calculation, one needs to use the conjugation formula 

u l B (0)q~ Bo UabW) = q~ Bo | [ 1 - (1 + q~ 1 )9 2 q~ B ° +d 2 + < ? “ 1 0 2 g“ (Ao+So) 

Y - <A„ +B 0 )/2 Ji-Q 2 q -B oAB+ -e q - 1 (l-q)A + B_ Jl-9 2 q~ B ° q “ (A ° +So)/2 1 


-6(l-q)q- 

which can be obtained in a fashion similar to the way 0831 ) was derived, as well as the formula 


I > 


U\ B (0)q~ Bol2 B - U AB (9) = q~ Bo/2 B _ yj 1 - 9 2 q~ B ° -9q~ B °A-q ~ A ° 12 , 
and its complex conjugate, which are easily obtained from the q-BCH relations. The final result is 

(93) g^o+^o = g Ao+So |l_(l + g- 1 )0 2 g- Co + ( /) 2 J 

+ q _1 0 2 g A °“ Co [ 1 -(1 + q~ l )0 2 q~ B ° + 9 2 + < ? “ 1 0 2 g“ (Ao+So) 

- q _1 0 2 0(l - q)q (Ao ~ Bo ~ 2Co),2 \/l-e 2 q- B »A-B + - q^^Od - q)A + B-\J 1 - 92 q -B 0q (Ao-B 0 -2C 0 )/2 

- 0(1 - q)q (2Ao+B<,)/2 B 1 - d 2 q~ B °q~ C ° /2 \/1 - 0 V C °C+ 

-Co/2 


0(1 - qW 1 - d 2 q- B °B + q i2Ao+Bo),2 C- Jl-(p 2 q~ c »q 


+ 00q _1 (l- q)A-q A ° /2 q~ Ca ' 2 1- (f> 2 q~ c °C+ + 00q _1 (l- q)q A ° /2 A + C_y l-(p 2 q~ c °q~ c ° 12 . 
As a consequence, we have the following 7-point eigenvalue equation 


(94) g»i+» 2 S^ B2 (*i,* 2 )= [gM(l-(l + g- 1 )0 2 q 


X1+X2-N + <p2 


+ q~ 1 (p 2 q 2xi+X2 ~ N 


(l - (1 + q” 1 )*) V* 2 +0 2 + q- 1 e 2 q-^ +x ^] ] E<*J, 2 (*i,* 2 ) 


- q~ 2 d(f) 2 q ( ' Sxi+Xi ~ 2N ^ 2 \/(1 - q X2+1 )(l - q x i)(l - 02 g -(**+l)) H^„ 2 (*i - l,x 2 + D 


- q- 1 00V 3X1+X2 - 2JV)/2 J( 1 - q X2 )(l - q x i +1 )(l - 0 V X2 ) S W ( Xl + 1,* 2 - 1) 


- 0 ? 

-0q 


(3*i+2* 2 -iV-2)/2 


(3xi+2x2-AT)/2 


(1 - q X2 )(l - 9 2 q~ x A( 1 - q^- x i- x 2+l)(i _ 02 g * 1+ * 2 -iV-i) EgJ, 2 (xi,x 2 - D 
(1 - q X2+1 )(l - 02 g -(* 2 +D)(i _ g JV-x 1 -x 2 )(i - (p^qXi+xi-N) E^* re2 (xi,x 2 + 1) 


+ q“ 1 00q (2xi+a:2 “ iV “ 1)/2 y / ( 1 - q x i)(l - - 02^+^-iV-l) S^J, a (*i - l,x 2 ) 


+ q _1 00q (2xi+X2_iV+1)/2 \/( 1 - q Xl+1 )(l - g iV-xi-*B)(i_02 g *i+xs-lV) E^fo + l,x 2 ). 
In terms of the polynomials, this relation amounts to 

(95) q rei+ ' l2 K„ li „ 2 (xi,x 2 )= [q Xl+X2 (l-(l + q- 1 )0 2 q Xl+X2 - iV + 0 2 ) 

’ (l - (1 + q~ 1 )9 2 q~ X2 +9 2 + q _1 0 2 q“ (xi+X2) ) 


+ q _1 0 2 q 2xi+X2 ^ ( 




+ 0"q x i- 2 "-- 1 )!- q Xl )(l -9 2 q~ X2 ~ 1 ) K rei)7l2 (xi - l,x 2 + 1) 


t2-*i+£2 - iV-l/ 

' J. — 1/ A-*- — u ^ — -*->• 

+ 9 2 <p 2 q 2xi ~ N1 (l - q X2 ) K nii „ 2 (jci + l,x 2 -1) 

+ q Xl (l - q X2 )(l - 0 2 q Xl+X2 “ iV - 1 ) K ni; „ 2 (xi,x 2 - 1) 

- 0V 1+X2 (1 - q Xl+X2 ~ N )(l - 9 2 q~ X2 ~ 1 ) K„ 1>B2 (*i ,* 2 + 1) 


+ (1 - q Xl )( 1 - 0 2 q Xl+X2 “ iV - 1 )- l,x 2 ) -0 2 0 2 q Xl - 1 (l - q x ^~ N ) K reij „ 2 (xi + l,x 2 ). 
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The eigenvalue equations < f88l > and 1951 ) coincide with those of Ifl6ll . 


5.2 Recurrence relations 

We shall now derive recurrence relations for the matrix elements S^,re 2 (xi,x 2 ) as well as for the 
two-variable polynomials K„ 1>re2 (xi,x 2 ). 


5.2.1 First recurrence relation 

To obtain a first recurrence relation for the matrix elements ’E^ ri 2 {x\,X 2 ), one can start from the 
eigenvalue equation ( 1871 ) and apply the duality relation ( l29l ) to exchange the roles played by the 
degree indices n\, n 2 and the variable indices x\, x 2 . Applying d29l ) on d87l ) and performing the 
change of variables and parameters 

X 2 ^n, 2 , N-n\ -n 2 <-*• N-x\ —X 2 , 6 *—(p, 

one finds 


?“* 1 “** E S 5 » a (*i.* 2 ) - -H~ (ni+2n2 ~ N)/2 \/( 1 - qW-'u-« 2 + 1)(1 _ g» 2 )(i -E^ ) ra2 _ 1 (xi,x 2 ) 

_j_ ^—(^i+2?i2+l) 


q n2+1 ( 1 + <P 2 ) + (f> 2 q ni+n2 ~ N - (1 + q) <p 2 \ S^J, 2 (xi,x 2 ) 


_ ( j )q -(n 1+ 2n 2+ N-2)/2^J il _ q n 2+ l )(1 _ q N- ni -n 2){1 _ 02 g - B2 -l) -^^(xi^). 

In terms of the polynomials, one finds 


(96) q“ Xl “ X2 K ; 


n\,n 2 


(xi,x 2 ) = q 


_ —(n.i+2/ii + l) 


? re2+1 (l + (p 2 ) + cp 2 q n 1+n2 ~ N - (1 + q) (p 2 


K 


ni,n 2 


(xi,x 2 ) 


+ q 


—(tzi + 1)/ 


(1-q 2 )(l-g 


ni+n 2 -N-l\ 


)(1 -cp 2 q re2 )K rei; „ 2 _i(xi,x 2 )-g 


— (Sn2~^fii + 1 ) l 2 , 


0 ^-ni,ri2+ l(Xl,X 2 ). 


5.2.2 Second recurrence relation 

To obtain the second recurrence relation, one applies the same procedure as above on the eigenvalue 
equation dMl ) with the result 


(97) q“ Xl S^ > ) „ 2 (xi,x 2 ) = q- (2ni+n2+N+1)l2 e<p^fa - q^)(l - g JV-m-n2+i) (1 - Q 2 q~^) S^ 1>n2 (xi,x 2 ) 

- q- &ni+naV2 0y/( 1 - q n * +1 )( 1 - q n ')a-e 2 q- n i)a-(p 2 q- n *- 1 ) S^ 1>Bg+1 (xi,x 2 ) 

- q- ll2 q- (2n ^ n ^Q 2 (P\j (1 - q-^)( 1 - q ni+n 2 -N-l )(1 _ ^2 q -n 2) E^ 2 l (xi,x 2 ) 


+9 


-2rei 


g rei (l+0 2 -q“ 1 ' ni (l + g)0 2 )+q' 1 ' re2 0 2 (l+0 2 + g rei ' iV 'V 2 -9“ n2 “ 1 (l + 9)0 2 ) 


“ re T,™2 (x l> X 2) 


—(3aii+ 4/12+-A^+4)/2/i2 


_ q , • -- 2 0^c/>\/ (1 - q^-«l-^ 2 )(l - g «2 + l)(l - (pZq-n. 2 - 1 ) H ( ^ j ) re2+1 (xi,x 2 ) 

_ q -(3n 1+ n 2+ l)/2 e ^ a _ g „ a)(1 _ q n 1+ l ){1 _ g2 q -n 1 -l )(1 _ -^^(x^X,) 


+ q (2ni+n2+N+3)/2 6(PyJ( 1 - q' 1 ^ 1 ))! - g tf-*i-» a) (i - dV^" 1 ) 1>ng (*i,x 2 ). 
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In terms of the polynomials, one finds 


(98) g-* 1 K B1 , B2 (* 1 ,* 2 ) = <?“ 2rei {g rei (l + e 2 - q~ n ^\l + q) 6 2 ) 

+ q-^d 2 (1 + (P 2 + q ni - N - 1 (p 2 - q-'^Hl + q)cp 2 ) j K Bl>B2 (xi, * 2 ) 

+ q~^ +2) (p 2 (1 - g rei )( 1 - g re i+™2-iV-l) (1 _ 0 2 g -«i) K ni _ 1)B2 (xi,x 2 ) 
+ q~ 2 ^ + 1 ) (p 2 (1 - g rei )( 1 - 0 2 9 - Bl )K Bl _i >Ba+ i(*i,* 2 ) 

+ g“ (2rei+2) 0 2 d - g re2 )(l - q ni+ri2 ~ N ~ 1 )(l - <p 2 q~ n *)K nun 2 _ 1 (x 1 ,x 2 ) 
- q-^ n 2+ 2 n i+ 2 '>g 2 c p 2 K„ 1;B2+ i(xi, x 2 )- g“ (3,ll+n2+1) 0 2 K Bi+ i iB2 (xi,x 2 ) 

+ 0 2 0- 2 g re 2-3rei-2 (1 _ g n 2)(1 _ 0 V" 2 ) 

Remark. Let us note that the structure relations with respect to the variable indices xi,x 2 derived 
in Subsections (4.3) and (4.4) cannot be obtained by the same duality argument from the structure 
relations with respect to the degree indices ui,u 2 given in Subsections (4.1) and (4.2). 


6 Conclusion 


In this paper, we have constructed an algebraic model for the two-variable quantum q-Krawtchouk 
polynomials. We have shown that these polynomials arise as matrix elements of a three-dimensional 
q -rotation operator expressed in terms of q -exponentials that act on the states of three independent 
q -oscillators. We have explained how this model can be used to derive the main properties of the 
polynomials: orthogonality relation, duality property, structure relations, eigenvalue equations and 
recurrence relations. 

It is apparent that the approach presented here can be extended to any number variables; the 
only essential difficulty lies in the complicated notation required to describe multivariate polynomi¬ 
als. Nevertheless, let us exhibit how the multivariate quantum g-Krawtchouk polynomials of Gasper 
and Rahman arise in the present framework. Consider d + 1 independent copies of the q -oscillator 
algebra with generators A®, A 1 ,' 1 and commutation relations 

[A^A^'W^A^, [A ( l ) ,Af] = S lJ q A o, i = l,2,...,d + l. 

Let {ni,n 2 .. .,nd\ and {xi,x 2 ,.. .,Xd) be two sets of non-negative integers such that Lf =1 n-d^N and 
X! k-i x k ^ N, where N is a positive integer. We define n = (ni,n 2 ,.. .,na+i) with rid+i = N - Lf =1 n-d, 
and similarly for x. We shall use the notation: 

ly*l = yi+y2 + ---yk- 


Note that by definition, we have ly^+il = N and |yol = 0. Consider the representation of the algebra 
generated by the d + 1 independent q -oscillators on the states 

d+1 

|n)jv = 01 n k ), 

k=i 


and defined by the following actions of the generators on the factors of the direct product: 


(99) Af\m) = 


‘l-qUi+i 

—T -l»i + l>, A_! |n,;> = 

1-g 


1 - q Ui a) 

— -|n;-l>, Ayim) = m\ni). 

1 -q u 
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We introduce the q-rotation operator in the (i,j) plane which is defined as 


U lJ (9) = ef [e 2 q- A ' 


(i) 

0 \e t 


(6a-q)q- (A o +A o >)/2 AfA^} 

x E q (-0(1 - g)g _(A o ,+A o >)/2 A ( _ i) AV’ ) ) E^ [~9 2 q~ A o’). 


Let H n (x;0i, ■■■ ,6d) be the following matrix elements: 


H n (x;0i,-” ,0 d ) = Ar(n| nf=o Ud-k,d-k+i(8d-k) \*)n- 

A direct calculation shows that in terms of the one-variable matrix elements d 10b . the matrix ele¬ 
ments H n (x;0i,--- ,9d) have the expression 


d 

H n (x;0i,• • • ,6d) = [I £ 


k=i 


( I x* + 1 1 — | n* -1 1 ) 
ra*,|x*|—|n*_i| 


(6k). 


In view of dill ), this means that these matrix elements are proportional to the multivariate quantum 
q-Krawtchouk polynomials K n (x;0i,.. .,9d',N) which read 

d I i 

K n (x;0i,...,0 d ;AO= Yl k nA \*k I - l n *-il, ^, l x A:+il - l n *-il; O' , 
k=l l u k 

where k n {x,p,N\q ) are the one-variable quantum q-Krawtchouk polynomials given by ©. With the 
above expression for Z n (x;0i, • • • ,9d), one can proceed with the determination of the properties of the 
multivariate q-Krawtchouk polynomials. 

As explained previously, the results presented here can be viewed as a q-generalization of IflOH (in 
the special Tratnik case), where an algebraic interpretation of Griffiths’ multivariate Krawtchouk 
polynomials was given in terms of the rotation group acting on oscillator states. The approach of 
IfTQH was generalized in |6j and 0, where algebraic interpretations of the multivariate Charlier 
and Meixner polynomials involving the Euclidean and pseudo-rotation groups were found. In light 
of these results, it would be of interest to generalize the approach developed here to also obtain 
algebraic models for the multivariate q-Charlier and q-Meixner polynomials. 
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